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The anisotropic motion of an interface driven by its intrinsic curvature or by an external field is 
investigated in the context of the kinetic Ising model in both two and three dimensions. We derive in 
two dimensions (2d) a continuum evolution equation for the density of kinks by a time-dependent and 
nonlocal mapping to the asymmetric exclusion process. Whereas kinks execute random walks biased 
by the external field and pile up vertically on the physical 2d lattice, then execute hard-core biased 
random walks on a transformed Id lattice. Their density obeys a nonlinear diffusion equation which 
can be transformed into the standard expression for the interface velocity v = M[(y + 7") ft + H], 
where M, 7 + 7", and ft are the interface mobility, stiffness, and curvature, respectively. In 3d, 
we obtain the velocity of a curved interface near the (100) orientation from an analysis of the self- 
similar evolution of 2d shrinking terraces. We show that this velocity is consistent with the one 
predicted from the 3d tensorial generalization of the law for anisotropic curvature-driven motion. 
In this generalization, both the interface stiffness tensor and the curvature tensor are singular at 
the (100) orientation. However, their product, which determines the interface velocity, is smooth. 
In addition, we illustrate how this kink-based kinetic description provides a useful framework for 
studying more complex situations by modeling the effect of immobile dilute impurities. 



I. INTRODUCTION 

Many bulk properties of polycrystals are strongly in- 
fluenced by the underlying microstructure. Much effort 
goes into predicting the motion of grain boundaries in 
response to a variety of driving forces. Depending on 
the nature of the grains, their boundaries migrate in re- 
sponse to applied stresses p] or magnetic fields Q , inter- 
nal forces associated with grain boundary curvature 
concentration gradients Q, etc. Successful models of mi- 
crostructure evolution must be supplied with the details 
of the ways in which the grain boundaries respond to the 
driving forces. 

Based on general conclusions of non-equilibrium statis- 
tical mechanics, one would expect the interface to have 
a unique mobility, i.e., a unique response coefficient to 
disparate driving forces. This conclusion was recently 
called into question by both experiments in polycrystals 
[H, |(| and simulations of Ising interfaces . These works 
observed drastically different shapes of shrinking grains. 
Grains shrinking under the influence of capillarity alone 
were roughly circular whereas the presence of other driv- 
ing forces resulted in strongly anisotropic shapes. This 
observation was most simply interpreted in terms of dif- 
ferent interfacial mobilities for different driving forces. A 
resolution of this apparent paradox in the Ising model 
0, which does not require a non- unique mobility, rests 
with identifying the crucial role of anisotropy in the cal- 
culation of the capillary driving force. This driving force 
is the strongly anisotropic interfacial stiffness i.e., the 
sum 7+7" of the excess free-energy of the interface 7 and 
its second derivative with respect to inclination, rather 
than 7 itself which is much less anisotropic. It turns 
out that the reduced mobility, i.e., the product of the 
capillary driving force and the bare mobility, is roughly 
isotropic, and therefore the grain shape is isotropic as 
well. Here we shed further light on the microscopic mech- 
anism for cancellation of the anisotropies of the interfa- 



cial stiffness and the interfacial mobility in both 2d and 
3d. 

The precise microscopic mechanisms responsible for 
the migration of grain boundaries are complex. However, 
there is hope that generic features near equilibrium are 
shared by a large class of models of moving interfaces. It 
is with this hope in mind we use the kinetic Ising model 
(KIM), introduced in Ref. @, as a proxy for studying 
grain boundaries. The KIM is defined by a collection of 
spins Si = ±1 on a lattice, a total energy which is a func- 
tion of this collection, and rules for dynamic evolution of 
the spins at some temperature (3 = 1/ kT. The energy in 
the presence of a magnetic field H is 

E=-jY,SiS j -H^s i , (1) 

(ij) * 

where the sum in the first term in Eq. Q is over pairs 
of nearest neighbors. Glauber dynamics [|| is one pos- 
sible scheme for evolving the collection of spins in such 
a way as to obtain correct distributions in equilibrium. 
This model is perhaps the simplest representation of non- 
equilibrium dynamics of interfaces. It can be used to 
explore effects of lattice anisotropy on the motion of do- 
main walls driven by magnetic field or capillary forces. In 
addition, domain nucleation and late stages of phase sep- 
aration can be addressed within the KIM. With simple 
modifications the KIM can be used to study the phe- 
nomenology of interface motion in the presence of mobile 
or quenched impurities. 

Much is known about the equilibrium behavior of Ising 
interfaces. For example, an exact expression for the in- 
terfacial free energy has been derived in 2d on a square 
lattice |0, Approximate expressions for this free 

energy and critical amplitudes in 3d have also been de- 
rived [121 IT^ I . The non-equilibrium behavior of the KIM 
is more complicated. Whereas several approximate ana- 
lytic results exist for the mobility of a domain wall in 2d 
Hi Ha IIS, little progress has been achieved is 3d. 
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Here we construct a simple and intuitive kinetic de- 
scription of low temperature domain walls in the KIM 
based on the kink degrees of freedom. The density of 
kinks is shown to obey a non-linear diffusion equation, 
which is equivalent to the law of anisotropic interface 
motion driven by curvature and/or an external field de- 
rived from the interface free-energy and mobility. It is 
important to emphasize that, in our kink description, we 
obtain the law of interface motion directly in the con- 
tinuum limit without these expressions as input into our 
calculation. Hence, our kink-based theory can be viewed 
as a direct microscopic derivation of the law of interface 
motion in the low temperature limit of the KIM, free of 
extraneous assumptions. Moreover, the kink-based de- 
scription is useful for analyzing more complicated situ- 
ations. We illustrate this point both by extending the 
analysis of anisotropic interface motion to 3d and by ex- 
amining impurity effects in 2d. 

Section[H]of this paper is devoted to the 2d KIM while 
the following Sec. lIIII extends our results to 3d. In section 
III Al we review the existing results concerning the KIM in 
2d with the focus on the non-equilibrium response of an 
interface to curvature and magnetic field. In Sec. Ill Bl we 
rederive the velocity of a curved Ising domain wall driven 
by capillary forces using kinks as the degrees of freedom 
responsible for the motion of the interface. This descrip- 
tion is accurate at low temperatures when the rate of 
nucleation of kink-antikink pairs is small. We obtain the 
shape of a shrinking Ising grain analytically in Sec. Ill CI 
In the following section lll Dl we illustrate the usefulness of 
the kink description by considering the influence of impu- 
rities on the grain boundary motion. In Sec. lIIII we study 
curvature driven motion in 3d near a high symmetry sin- 
gular orientation where the interface can be represented 
by a collection of terraces composed of kinks. This allows 
us to use the 2d analytic results to calculate the interface 
velocity and therefore the mobility tensor near this sym- 
metry direction. Finally, conclusions are given in section 
IV. 



II. TWO-DIMENSIONAL KINETIC ISING 
MODEL 



A. Low temperature expansion of the interface free 
energy and mobility 

Let us summarize the analytical results obtained so far 
for the KIM, focusing on the expressions that have been 
derived for the interface mob ility . In 2d, the exact in- 
terfacial free energy is known hOlhJ. For our purposes 
it suffices to write down the first two terms in the tem- 
perature expansion (enthalpic and entropic respectively). 
When the spins are arranged on a square (denoted by a 
□) lattice of unit lattice spacing, this energy is 

7 2 D d (</>) = 2J{c+s) + \ [c logc+s log s-(c+s) log(c+s)], 



where c = | cos 0| , s — |sin0|, and (f> is the inclination 
defined as the angle of the interface normal with respect 
to the (10) axis of the underlying lattice. 

When spins flip according to non-conserved Glauber 
dynamics |9(, the interface moves to minimize the free 
energy of the system which consists of the bulk and the 
interface contributions. Spohn [LH has derived the sharp 
interface continuum description of a domain wall in KIM. 
It follows from his derivation that the normal velocity of 
the interface v is the product of the mobility M and a 
driving force. From the continuum description it also 
follows that in the absence of magnetic field, the driving 
force is the product of the mean curvature of the inter- 
face k and the interface stiffness. In 2d, the stiffness is 
7 + 7", where 7" denotes the second derivative of 7 with 
respect to <j). Using a Green-Kubo perturbative formal- 
ism, Spohn obtained the interface mobility in the limit 
of small temperature and small driving magnetic field. 
The same result was obtained earlier by Barma [l5| us- 
ing a mapping of the dynamics of the low-temperature 
Ising interface to the one-dimensional exclusion process. 
Rikvold and Kolesik obtained analytical expressions valid 
for large fields and temperatures 16] . The leading term 
in the temperature expansion of the mobility diverges like 
1/T 



ft 



I sin 2(j>\ 



2r I cos</>| + I sin< 
/V3 sin (cos 0- 



V3' 



2r 



v/3 



(3a) 
(3b) 



where □ refers to a square lattice and A to a triangular 
lattice. In addition, t is the intrinsic time scale of the 
Glauber dynamics which is the inverse frequency of the 
attempted spin flips. The triangular lattice formula is 
valid in the <f> € [0, 7r/6] domain and can be extended to 
the other angles via an appropriate symmetry transfor- 
mation. 

The above expressions for the interface energy and 
mobility can be combined to arrive at the continuum 
(or mean field) low temperature equation of motion of 
the interface driven by curvature k and magnetic field 
H <C 1//3. The normal velocity of the interface is 

v 2d (<j>) = M 2d [K(72 d + 72d) + H] = M| d k + M 2d H, (4) 

where the reduced mobility on the square lattice is (see 
Appendix 1X1 for the triangular lattice result) 



M* d (^)^M 2 c d (7 2d + 72 ' d ) = 



r (| cos 0| + I sin0|) 5 



(5) 



ft 



(2) 



Note that the reduced mobility is roughly isotropic 
whereas the bare mobility is strongly anisotropic. In 
addition, the reduced mobility does not diverge in the 
T — > limit whereas the bare mobility does. This 
happens because the contribution of the enthalpic term 
2 J(c + s) to the interfacial stiffness evaluates to zero. 
Therefore, only the entropic term (which is proportional 
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to T) contributes to the stiffness. Moreover, the con- 
tribution due to the entropic term diverges at the high 
symmetry orientations whereas the bare mobility van- 
ishes at those orientations in such a way that the product 
of the two quantities produces a finite non-zero reduced 
mobility. This behavior is responsible for the nearly cir- 
cular shape of a shrinking grain on a hexagonal lattice in 
Ref. 0. 

B. Direct calculation of M* via the dynamics of 
kinks 

The basic law of interface motion embodied in Eq. Q 
is usually derived using a thermodynamic approach 
where the interface free-energy and mobility are com- 
puted separately. This approach, even though general, 
lacks intuitive appeal. Furthermore, it is not simply 
extended to more complex situations. It is therefore 
worthwhile to develop an alternative method for deriv- 
ing Eq. (JEJ directly from a microscopic picture without 
the need to compute the interface free-energy and mo- 
bility as intermediate steps. We develop such a method 
based on a low temperature description of the interface 
in terms of kinks. This simple microscopic picture, and 
the results obtained for the velocity of a shrinking grain 
in 2d, provide the basis for the subsequent incorporation 
of impurities and the derivation of an expression for the 
interface velocity in 3d. 



o 


o 


o 


o 


o 


o 


o 


o 


o 


X 


X 


X 




o 


o 


o 


o 


o 


o 


X 


X 


X 


X 


X 


X 




o 


o 


o 


X 


X 


X 


X 


X 


X 


X 


X 


X 



FIG. 1: Schematic description of an interface between up 
(crosses) and down (circles) spin domains on a square lat- 
tice. Only corner (boxed) spins can flip at low temperature. 
A flip of a corner spin corresponds to moving the kink left or 
right. 



When the temperature is low, j3 J 3> 1, the only al- 
lowed spin flips are those that do not increase the total 
energy. Therefore kink-antikink pairs cannot nucleate at 
the interface. Barma 0] observed that the interface be- 
tween the Ising domains can then be represented by a 
staircase of kinks shown in Fig. ^ The kinetics of kinks 
reduces to an exclusion process (asymmetric in the pres- 
ence of magnetic field) flTL fl8| . Even though steady state 
properties of this process (corresponding to a flat field- 
driven interface of a fixed inclination) are well known, 
little progress has been made analytically to describe the 
evolution of a non-uniform kink distribution correspond- 
ing to a curved interface. 

Let us define the ensemble average density of kinks 



p(x, t) and derive its evolution equation, which is equiv- 
alent to Eq. |0J. We outline the derivation here and 
relegate the details to Appendix [X] We focus here on 
the curvature driven motion while the appendix includes 
the effect of the magnetic field. Unimpeded by its neigh- 
bors, each kink executes a random walk corresponding 
to purely diffusive motion. Many kinks can "pile-up" at 
the same site but cannot pass through each other. Via a 
transformation which inserts an extra lattice site between 
every pair of neighboring kinks (illustrated in Fig.[2J), we 
map the dynamics of kinks onto the problem of Id ran- 
dom walkers which cannot occupy the same lattice site. 
The density of walkers for this symmetric exclusion pro- 
cess obeys a simple diffusion equation This is true 
because when two walkers collide, their indices can be 
exchanged (i.e. their "identities" switched) without af- 
fecting their density. Thus a collision can be viewed as 
the tunneling of the kinks through each other without 
affecting each other. The density of hard core random 
walkers is insensitive to this identity-switching transfor- 
mation and must therefore satisfy a diffusion equation. 
When transformed back to the original coordinate system 
in which kinks can pile up, the equation for p{x,t) reads 
for the square lattice (see Appendix^] for the triangular 
lattice version) 

where subscripts denote differentiation, F = —p x / (1+/?) 2 
is the flux of kinks and p = — 1/(1 +p) is the kink "chem- 
ical potential." Eq. © is a nonlinear diffusion equation 
with the diffusivity 1/(1 + p) 2 which is a decreasing func- 
tion of density. This reduction results from the fact that 
when more than two kinks occupy the same site, some of 
these kinks are completely immobile. Since the density 
of kinks is defined for interfaces inclined with respect to 
the (10) orientation, Eq. © has to be supplemented by 
boundary conditions which piece together different tt/2 
sectors of the grain boundary. 




FIG. 2: Mapping of the kink dynamics onto the symmetric 
exclusion process corresponding to non-overlapping random 
walks on a Id lattice. Note that kinks that pile up vertically 
in the physical 2d lattice do not overlap in the transformed 
Id lattice. See Appendix 1X1 for details. 

Geometrically, the density of kinks is the local slope 
of the interface with respect to the low energy (10) ori- 
entation. It is therefore easy to show that Eq. Q with 
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H = and Eq. © are equivalent (see Appendix lA"|l . 
Thus we derived the equation of motion for the interface 
without assuming the applicability of the continuum de- 
scription of the interface. Even though, for clarity, we 
have restricted our discussion above to motion by curva- 
ture only, we derive in AppendixE]the evolution equation 
for the kink density for general motion by both curvature 
and an external field, and show that it is equivalent to 
Eq. |gj). 

Neglecting thermal excitation of kink-antikink pairs al- 
lowed us to construct an equation for a single density of 
kinks p{x,t). In general local densities of kinks p+ and 
anti-kinks p_ must be considered. Each density obeys the 
non-linear diffusion equation © augmented by a source 
term proportional to exp(— 2/3 J), due to the creation of 
kink-antikink pairs, and a sink term proportional to the 
product p+p- due to the annihilation of kinks by anti- 
kinks. The local slope of the interface with respect to the 
low energy orientation is given by the sum p+ + p_ of the 
kink and anti-kink densities. Once the details of this two- 
density approach are worked out, a formal temperature 
perturbation expansion becomes possible since either p + 
or p_ is exponentially small in low temperature limit. A 
small mobility for the high symmetry orientations, which 
are immobile at zeroth order, will be the most important 
effect at the next order in the temperature expansion. 



C. Evolution of an Ising grain 



vanishes for these orientations. Thus a grain cannot grow 
beyond its initial size. Any growth process has to include 
the nucleation of kink-antikink pairs which is explicitly 
ignored in our description. 

The amplitude of the positive magnetic field H deter- 
mines the shape of the shrinking grain. When (3H ^> 
j3H c = 1/R, the second term in Eq. Q dominates. Note 
that for large grains this crossover magnetic field vanishes 
like 1/R. The shape of the grain shrinking under these 
conditions, shown in Fig. [2t>, is strongly anisotropic. 

When the applied field is much smaller then the 
crossover field, the evolution is controlled by the more 
isotropic reduced mobility and thus the shape of a shrink- 
ing grain is close to a circle. Even when initially the dy- 
namics is controlled by the magnetic field, the crossover 
to curvature dominated dynamics will happen when the 
grain shrinks to a sufficiently small size. In this regime, 
the grain shrinks in a self-similar manner (see Fig.[3J;). 



1. Self-similar evolution of the shrinking grain 

Here we restrict ourselves to the square lattice while 
citing the results for the hexagonal lattice in Appendix 
IB1 We also set the time scale r = 1. 



Let us now use the equation of motion to describe 
the evolution of an "Ising grain," i.e., an island of down 
spins in a sea of up spins. The reduction of the interfacial 
free energy and spin alignment parallel to a positive mag- 
netic field are the the driving forces for the grain shrink- 
age while negative magnetic field favors grain growth. 
The numerical method of solving J2J is described in de- 
tail elsewhere Q. 




FIG. 3: Shape of an evolving Ising domain in 2d. (a) Final 
stationary shape of a domain in sufficiently strong negative 
magnetic field, (b) Shape of a shrinking domain in a large 
positive field, (c) Self-similar shape of a domain shrinking in 
absence of magnetic field. 

When a sufficiently large negative is applied, the grain 
grows until it reaches a stationary configuration deter- 
mined by its initial shape (see Fig. This happens 
because the velocity of the interface in the direction of 
the low energy planes is never outward since the mobility 



£(t) 



FIG. 4: A 7r/4 slice of the grain defined by the dashed lines 
that converge in the center of the grain. The thick line is 
the grain boundary. The kink density p(x,t) is defined on a 
shrinking domain of width £(t). 

To compute the shape of a shrinking grain we need 
to specify the region in which our kink description holds 
and fix the boundary conditions at the edges of this re- 
gion. Since we expect the grain to possess four mirror 
planes inclined at 0, ±7r/4 and ir/2 with respect to the 
(10) plane, we will restrict ourselves to a 7r/4 wedge (See 
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Fig. UJ. The slope at the left edge of the wedge (x = 0) 
is p = and the slope at the right edge (x = £(t)) is 
p = 1 due to mirror symmetry around the tt/4 plane and 
the smoothness of the grain shape. Thus we are to solve 
Eq. JBJ) subject to the boundary conditions 

p(0,*)=0, p(£(t),t) = l. (7) 

The final ingredient in determining the grain shape 
is the shrinking rate. The slice width £{t) shrinks as 
the kinks at its right edge flow to the left with a flux 
F(£) = —p x (£)/A. Every time the kinks move two sites 
to the left, the width of the slice is reduced by 1 (see 
Fig. [3] for a visual explanation) , and therefore 



l(t) = 



F(£) _ Px (£) 



(8) 



-new 



^old_ 



FIG. 5: When the rightmost kink moves two steps to the left, 
the width of the domain is reduced by 1. 

As we mentioned above, in the absence of a magnetic 
field, the shape of the shrinking grain is self-similar. To 
prove this we seek a solution to the moving boundary 
problem defined by equations |JB}, J7J and (JSJ, which de- 
pends on space and time only through a combination 
C = x/£(t). Substituting this Ansatz into the expression 
for the shrinking rate (|HJ we obtain 



££ = - 



(9) 



where prime denotes differentiation with respect to £. 
Thus the the rate of change of the grain area A ~ £ 2 
under this self-similar evolution is constant as expected. 
The kink diffusion equation 10 becomes 



B(p' = 



(1 + P) : 



(10) 



with p(0) = 0, p(l) = 1, and B = -££ = p'(l)/8. The 
constant B 0.331491 is determined self-consistently by 
a shooting procedure. Fig.Elshows the comparison of the 
solution of (|10l) to the ensemble averaged Monte Carlo 
simulation of diffusing hard-core kinks with boundary 
conditions appropriate to the shrinking grain scenario. 
Quantitative agreement of the sharp interface result (0J 
with the Monte Carlo simulation of the full KIM was 
found in Ref. 0. 
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FIG. 6: Density of kinks in units of the inverse lattice constant 
plotted against the dimensionless scaled distance C, across the 
arc. The exact density for a self-similar shrinking grain ob- 
tained by solving 1101 is compared to that obtained via a 
Monte Carlo simulation of diffusing impermeable kinks. The 
dashed line shows for comparison the density of kinks in a 
circular arc on a square lattice. 

It is useful to recast Eq. Ill Oil in terms of the po- 
lar parameterization of the self-similar shrinking grain 
r((p,t) — \/2l{t) r2d(4>)- We chose to scale r((j), t) in such 
a way that r2d(7r/4) = 1. The kink density of a self- 
similarly shrinking grain is a function of <f> only 



r 2 d{(f>) sin0 - r' 2d ((f)) cos(f> 
r 2d (0) cos <p + r' 2d (4>) sin0' 



(11) 



The equation for the shape of the self-similarly shrinking 
grain in the polar coordinates can be integrated once to 
yield 



cos < 



r' 2d (<l>)(sm(f>- cos </>)] 2 



= rU4>) + 2r' 2 M-r'±M, (12) 

subject to r' 2d {0) — r' 2d {-K/A) = (by symmetry) and 
r2d(7r/4) = 1. One of these conditions is automatically 
satisfied for the value of B found above. 

Let us finally mention another analytic result concern- 
ing the grain shrinking rate dA/dt 



_dA 
~dt 



dcf>M* d (<f>) 



4, 



square, 



3-\/3, hexagonal. 



(13) 



These formulas (the square lattice result first appeared in 
Ref. [l9| ) are a simple consequence of the fact that only 
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the corner spins are allowed to flip (see Fig. [JJ. When 
a spin in a concave corner flips, the area of the grain 
increases by 1 (square lattice). And vice versa, when 
a kink in convex corner flips, the area is reduced by 1. 
Since the probabilities of all allowed spin flips are the 
same and the number of convex kinks on a square lattice 
is greater than the number of concave kinks by 4 (due to 
Hopf's theorem which states that the rotation index of 
a simple curve is 1), we arrive at lpH!|) . Ref. [8j checked 
that the shrinking rate on the hexagonal lattice is indeed 
3V3 « 5.196. 




FIG. 7: The number of the outside corner spins (diagonal 
Crosshatch) for any domain is always smaller by four than 
the number of the inside corner spins (horizontal-vertical 
Crosshatch). The shrinking rate of this domain, computed 
from the difference between the inside and outside corner 
spins, is thus known exactly in the low temperature limit. 



D. Drag by immobile impurities 

The kink picture of the low-temperature grain bound- 
ary dynamics is useful in understanding the effect of di- 
lute immobile impurities. We model the interaction of 
the grain boundary with interstitial impurities by defin- 
ing a variable 9 mn on the dual lattice sites. 9 = 1 when 
an impurity is present and otherwise. The impurities 
are randomly positioned on the dual lattice and do not 
move. The interaction of the impurities with the spins is 
introduced via an additional term in the energy 

m,n 

where S( mn } is the total magnetization of the Ising spins 
nearest to the impurity located at site (m, n) of the dual 



lattice and the sum is over all the dual lattice sites. 
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FIG. 8: Illustration of the energy landscape resulting from 
the interaction of a kink with a single fixed impurity (large 
solid circle) on the dual lattice. The up (down) spin domain is 
above (below) the interface represented by a thick black line 
along the dual lattice. The vertical segments denote three 
positions of the kink. The energy (in arbitrary units) corre- 
sponding to these three position is shown schematically below. 
When the impurity is positioned such that the top edge of the 
kink passes the impurity from left to right, the total energy of 
the system decreases (a). The opposite is true if the impurity 
is on the lower edge of the kink (b) . 

Figure [S] explains graphically that, depending on its 
position, an isolated impurity provides either a left or a 
right directed short range force acting on a kink. In ad- 
dition to this force there is a two-kink effect which makes 
kink pile-ups energetically favorable when they occur on 
the impurity site. 

So far we considered only positively charged impuri- 
ties 6 = +1. The interface is attracted to these impu- 
rities. Negatively charged impurities with 9 = — 1 repel 
the interface. However, the qualitative picture of the 
kink-impurity interaction presented in Fig. [S] still holds. 
The only difference is that the effect of the negative im- 
purity on the top edge of the kink is equivalent to the 
effect of the positive impurity on the bottom edge and 
vice versa. In the limit of high density of impurities, 
additional effects due to the interplay of positively and 
negatively charged impurities become important. For ex- 
ample, a row of alternating positive and negative impu- 
rities perpendicular to the interface pulls the interface 
along in one direction or another due to the ease of nu- 
cleating kink-antikink pairs. Additional phenomena arise 
when both positively and negatively charged impurities 
are presented. Exploring these these phenomena is out- 
side the scope of this article. 

The diffusing kink picture is especially simple when the 
impurities are dilute and £> J> kT. In this limit the 
kinks diffuse only downslope in the static energy land- 
scape produced by then impurities. When impurities are 
dilute, this energy landscape consists of a number of flat 
terraces bound by steep vertical cliffs or walls. The kinks 
diffuse and fall down cliffs until they fall onto a terrace 
which is bound by walls on both sides. The kinks become 
trapped on this terrace. In the long time limit, the den- 
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sity of the trapped kinks becomes uniform on this terrace. 
This means that the slope of the piece of grain boundary 
which corresponds to this terrace is a constant given by 
the density of the trapped kinks. This density depends 
on the initial distribution of kinks and impurities and can 
be anything. Therefore, in this limit, the grain boundary 
is pinned and consists of a series of flat facets of random 
length and inclination. 

Figure El presents results of the Monte-Carlo simula- 
tion of the low temperature 2d Ising model in the pres- 
ence of strong dilute positively charged immobile impuri- 
ties. What is shown is the time-averaged location of the 
boundary between the spin up and spin down domains. 
The boundary is pinned and consists of straight pieces 
of random length and orientation. Impurities are located 
at either end of each such facet. This result supports our 
qualitative picture. 

Addition of magnetic field introduces yet another en- 
ergy scale H into the picture. When H -C e, the bound- 
ary is pinned. The shape of the pinned facets depends on 
the relative size of 1 / (3H and length of the pinned facet 
L. In equilibrium, the kink drift due to the magnetic field 
is balanced by the diffusion due to curvature. Thus small 
facets for which L -C l/f3H, remain straight. Conversely, 
when L ^> 1/{3H, the long facets look like the corners of 
a droplet expanding in co-aligned magnetic field whose 
shape is given in Fig. EJJi. 

Strong positive impurities (e J) in the bulk of a 
domain of aligned spins always have two spins near then 
that are anti-aligned with the rest of the spins in that do- 
main. Thus impurities serve as nuclei for the formation 
of droplets of the phase of spins aligned with the applied 
magnetic field. Conversely, strong negatively charged im- 
purities favor alignment of the nearby spins and thus can 
inhibit nucleation of the phase favored by the application 
of magnetic field. 




FIG. 9: Time averaged shape of a grain boundary pinned 
by strong e — 5 impurities which occupy 1% of the dual lat- 
tice sites is a collection of straight segments obtained via a 
Monte-Carlo simulation of the 2d Ising model with an addi- 
tion energy given in Eq. 11411 . 



III. THREE DIMENSIONAL ISING MODEL 

Whereas a curve on a plane can be characterized by a 
single scalar curvature, a smooth surface embedded in a 
three dimensional space is characterized by a rank two 
tensor L a p ({a, (3} — 1,2). This tensor is called the 
second fundamental form or the Weingarten map or just 
the curvature tensor. The trace of this tensor is the mean 
curvature, while its determinant is the Gaussian curva- 
ture. This tensor is defined at some point P by selecting 
an orthogonal coordinate system x a in the tangent plane 
at P and writing 

dh 

L a !3 = t a ■ — , (15) 

where t a are the unit tangent vectors, and h is the unit 
normal to the surface. 

The reduced mobility of a two-dimensional interface 
is also a rank two tensor M* a which when contracted 

ap 

with the curvature tensor yields the normal velocity of 
the interface 

u,f3 

The reduced mobility tensor depends on the scalar bare 
mobility M%d (found, for example, by measuring the 
speed of a driven flat interface) and the interfacial free 
energy 7 and its derivatives. In the neighborhood of the 
point P the normal h is specified by its the deviations 
ipi and if2 from the normal at P in the directions t\ and 
i 2 - The free energy is a function of the normal h and 
therefore, in the neighborhood of P, a function of these 
angles 73d ("^l, y^)- The reduced mobility tensor is then 
defined as 

M ^M 7 ^ + s£fc-) (i7) 

Since our KMC simulations show that the shape of a 
3d shrinking grain is even closer to a sphere than a 2d 
shape to a circle, this reduced mobility tensor is nearly 
isotropic. This isotropy allows us to predict the 3d grain 
shrinking rate (defined as the rate of change of the 2/3 
power of its volume S = V" 2 / 3 ) by calculating the ve- 
locity u 10 o of the shrinking grain boundary at the (100) 
orientation. We will first estimate this velocity within the 
terrace-step-kink description of the vicinal surface. We 
then derive an exact expression for this velocity within 
the continuum limit. 



A. Shrinking terrace view of the dynamics near 
(100) plane 

The low temperature interface can be described within 
the the terrace-step- kink (TSK) model [2(j. When steps 
are far apart, each step obeys the dynamics of a 2d grain. 
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FIG. 10: Grain shape near the 100 plane can be viewed as a 
collection of terraces 



If the 3d grain is a sphere of radius R, it is described near 
the (100) orientation by a stack of circular terraces of in- 
creasing radii ri(i), rzit), etc (See Fig. llO|) . Because these 
steps are part of the spherical grain, they are related via 



R 2 



i = Jr 2 -' 



(18) 



Solving this equation for R 2 , differentiating with respect 
to time and using the exact expression for the 2d grain 
shrinking rate (|13f) rifi = T%r% = — 2/7T we obtain 
RR = — 2/tt. Therefore, within the spherical grain ap- 
proximation, the grain shrinking rate is 



S 



4n \ 2/3 4 



-3.309, 



(19) 



which is in reasonable good quantitative agreement with 
the shrinking rate found by KMC simulations S — 
-3.335 ±0.001. 



B. Low temperature expansion of interface free 
energy and mobility near the (100) plane 

Little analytical progress in deriving equilibrium and 
kinetic properties of the 3d KIM has been achieved to 
date. A mean field expression for the free energy of the 
TSK model neglecting step-step interaction was obtained 
by Gruber and Mullins Holzer and Wortis H2| cal- 
culate the free energy near the (100) plane in the more 
controlled diagrammatic temperature expansion. It is 
sufficient for our purposes to keep only the leading term 
in the expansion 



73d (Q,(t>) ~ 72d(0), 



(20) 



(|20|l is simply a statement that the free energy of a vic- 
inal surface is composed of the free energies of the steps 
(considered to be non- interacting) . 

In the same non-interacting step approximation, the 
bare mobility of the vicinal surface is 



M 3 d(0,0)«0Af2d(<£). 



(21) 



C. Speed of the (100) orientation of the shrinking 
grain 

The expressions for the interfacial free energy and the 
mobility in the vicinity of the (100) plane allow us to 
calculate the reduced mobility at this orientation as well 
as the grain shape in the neighborhood of the point P 
where the normal is in the z direction. 

Let the shape of the grain in cylindrical coordinates be 
z(r,4>) — r 2 1 r| d (4>) . The shape of the z = 1 section of the 
3d grain is r — r2d{4>) in polar coordinates. A circular 
terrace, i.e., r2d(4>) — V^R, corresponds to a sphere of 
radius R. We choose this suggestive parameterization of 
the 3d shape with the foresight that r2d(<f>) will turn out 
to be identical to the shape of the self-similarly shrinking 
2d grain. This is not surprising in view of the shrinking 
terrace picture of the 3d grain evolution near the (100) 
orientation of the previous subsection. 

Given the shape of the 3d grain, we can compute the 
normal of the surface n(r, <j)) and the curvature tensor 
L a p(r, 4>). The reduced mobility tensor M^p(h), in turn, 
can be calculated using the free energy expression i|20|) 
following the prescription (|17|) . Their contraction is the 
local velocity of the interface v(r, cf>). We will not write 
down the full expression for v due to its uncnlightcning 
complexity. It's r — > limit has to be independent of the 
direction of approach </>. We obtain 



lim ii(r, 4>) 



2 [rlM + 2r 2d (0) 



-r 2d (0)(sin0-cos <f)] 
= const, (22) 



subject to the smoothness constraints r' 2d (0) 



'2d 



(0) 



0. The function r2d{4>) which satisfies the above equation 
and constraints is precisely the self-similar shape of a 2d 
shrinking grain. If we set the size of the 3d grain by 
choosing r2d(7r/4) = V^R, we arrive at v = —2B/R and 
hence (again assuming a spherical shape to estimate the 
volume) 



2/3 



AB : 



-3.445. 



(23) 



where is the angle between the normal to the inter- 
face h and the z-axis (which is the normal to the (100) 
plane) and <f> is the angle between the projection of h onto 
the x-y plane and the x-axis counted clockwise. Formula 



Since the diameter of the self-similar shape at (j> = is 
slightly smaller than its diameter at <fi — 7r/4, a better 
approximation would have been to scale r2d in such a 
way that at the intermediate angle r2d(7r/8) = \/2R. 



9 



IV. CONCLUSIONS 

In summary, using a kink-based description, we have 
derived directly from a microscopic model (low tem- 
perature KIM) a continuum evolution equation for the 
anisotropic motion of a simple interface, and we have 
shown its equivalence to the standard phenomenological 
law of motion by curvature. We have illustrated with 
the example of dilute impurities that this kink-based ki- 
netic description provides a useful framework for study- 
ing more complex situations. By extending this descrip- 
tion to 3d, and by exploiting our 2d result for the self- 
similar dynamics of shrinking terraces, we have obtained 
the velocity of a curved interface near a singular orienta- 
tion. We have shown that even though the interface stiff- 
ness tensor and the curvature tensor are singular at the 
(100) orientation, their product, which determines the 
interface velocity, is smooth. Furthermore, this velocity 
is consistent with the one predicted from the 3d tensorial 
generalization of the law for anisotropic curvature-driven 
motion using known expressions for the interface free en- 
ergy and bare mobility. 

Our kink-based derivation of a continuum equation of 
interface motion highlights the microscopic mechanism 
for the remarkable isotropy of the reduced mobility in 
both 2d and 3d and thus the shape of grains shrinking un- 
der the influence of capillarity alone. The reduced mobil- 
ity is a product of the interfacial stiffness and the interfa- 
cial mobility both of which are strongly anisotropic. The 
isotropy of the reduced mobility is therefore a result of 
the cancellation of anisotropics of the interfacial stiffness 
and interfacial mobility. The microscopic reason for the 
cancellation is purely geometric in origin. The number 
of geometrically necessary kinks, and hence the configu- 
rational entropy of the interface, varies rapidly with in- 
clination near low-energy/low-mobility orientations, but 
slowly near high-energy/high-mobility interfaces, where 
the density of kinks is high. Since the leading order con- 
tribution to the interfacial stiffness comes from configu- 
rational entropy, stiffness is high where mobility is low 
and vice versa. The cancellation of anisotropics leads to 
roughly isotropic reduced mobility. Therefore the shape 
of a shrinking grain can appear isotropic or anisotropic 
depending on whether driving forces other than capillar- 
ity are present. The bare mobility of the interface is, 
however, independent of the nature of the driving force. 

An interesting prospect for the future is to extend 
this kink-based theoretical description of interface mo- 
tion to realistic, and more complex, grain boundaries 
where kinks have the character of dislocations. Work 
along this line is presently in progress. 

We thank Bernard Derrida, Mikhail Mendelev, An- 
thony Rollet, and David Srolovitz for valuable discus- 
sions. This research is supported by U.S. DOE through 
Grant No. DE-FG02-92ER45471 and funds from the 
Computational Materials Science Network. 



APPENDIX A: DERIVATION OF THE KINK 
EQUATION OF MOTION 

In general, kinks comprising the grain boundary are 
characterized by their width b, which is the distance 
of the closest approach of two neighboring kinks, their 
height d, and the length of the steps of their random 
walk a. For example, on square lattice, d = a, the lat- 
tice constant and 6 = 0, while on a triangular lattice 
d = av3/2 and b = a/2. In the continuum limit, we 
define the density of kinks p(x, t) and seek its evolution 
equation in some fixed domain x £ [x-i,, Xr]. Since neigh- 
boring domains contain anti-kinks, absorbing boundary 
conditions must be imposed p{x\ Jl t) = p(xp. , i) = 0. 

A special case of this problem a = b describes random 
walkers in Id which cannot occupy the same site. When 
a magnetic field is present, the random walk is biased and 
the problem can be map ped onto the well studied asym- 
metric exclusion process [jjl |2J| . To map the problem of 
finding the evolution of the kink density p(x, t) onto this 
special problem, we insert a space c = a—b between each 
pair of adjacent kinks, as illustrated for a square lattice 
in Fig. 13 The resulting kink density R(£, t) e [0, l/a] 
is defined in a different domain £ € [£l(£),£r(£)]- In the 
presence of magnetic field H, this kink density satisfies 
the equation [171 l24| (subscripts denote differentiation) 

R t = DR a + a[R(l - aR)}^ = - J s , (Al) 

with 

J(£, t) = -DR i - aR{\ - aR), (A2) 

where D = a 2 /2r, a = cl(3H/t, t is the Monte Carlo time 
step, and J is the flux of kinks in the moving ^-domain. 
Note that this equation is identical to Burger's equation 
after elimination of the drift term aR^ by transformation 
to a moving frame. As kinks annihilate at the boundaries 
of the x-domain, the ^-domain shrinks. Each kink that 
leaves the x-domain, decreases the ^-domain by c. This 
implies that the boundaries of the ^-domain move with a 
velocities proportional to the current of kinks our of the 
domain 

Sl = -cJ(&(*),*), £K = -cJ(&(t),t). (A3) 

The equations for the motion of boundaries l|A3fl . to- 
gether with the absorbing boundary conditions 

iZ(& (*),*) =R(£ R (t),t) = 0, (A4) 

completely define the problem of diffusing kinks in the 
£- domain. 

The mapping is inverted via 

£(x, t) — = x — xl + c / p(x',t)dx. (A5) 

J XL 

At some fixed time we can write 

R(£, t)d£ = p(x, t) dx, (A6) 
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since both expressions give the number of kinks in the 
same physical interval. Using (|A5() we obtain 

T>(t 4\ P^ X ^ I t\ R (^t) 



obtain 



1 + cp(x, t) ' 



(A7) 

This relationship (|A6|I allows us to invert (|A5|) to obtain 



-*l = -c / R(?,t)d?. (A8) 

It is now only a matter of carrying out the chain rule to- 
gether with the boundary conditions (|A4|I and the trans- 
formation HA5fl to obtain 



Rt = 



Px 



(l + cp)3- 
Pi 



D 



(l + cp) 2 (l + cp) 5 
p x cp(l - bp) 



—a ■ 



(i + c P y 



R 



3cpl 



(l + cp) 4 (l + cp) 



5 " 



(A9a) 

(A9b) 
(A9c) 



Thus we obtain the non-linear diffusion equation for 
p(x,t) which reads 



Pt = D 



where 



Pxa 



2cD- 



pi 



(1 + cp) 2 (1 + cpf 

p x (l - 2bp - bcp 2 ) 



a ■ 



-F x , (A10) 



t? n P x I ~bp 

F = - D T^T^- ap —p> (A11) 

is the flux of kinks in the fixed domain which vanishes at 
zero kink density. 

Using the relationship of the local slope and kink den- 
sity tan0 = h x (x,t) = dp(x,t) and the expressions for 
the normal interface velocity and curvature 



(1 + ^)1/2' 



(1 + ^)3/2 



(A12) 



we can compute the bare and the reduced mobilities from 
the normal velocity of the interface v n = M* k+MH . We 



M* = 



D {3 A sin <j> (cos — p sin <j>) 

(cos0 + sin 0)2 ' t cos </> + ^ sin ^ 

(A13) 

where v = c/d, p = b/d and A = a/d are geometric 
factors. These expressions are valid for (j> £ [0, 7r/4] for 
the square lattice and for € [0, tt/6] for the triangular 
lattice. 

APPENDIX B: SELF-SIMILAR SHRINKING 
GRAIN ON A HEXAGONAL LATTICE 

The symmetry of the hexagonal lattice allows us to 
solve for the shape of the self-similarly shrinking grain in 
a 7r/6 wedge. We present yet another way of obtaining 
this shape. Let the points on the boundary be labeled 
by (f>, the azimuthal angle 4> £ [0,7r/6]. Let 9(4>) be the 
local slope and r 2c i(0) t ne radial distance from the center 
of the grain. The shrinking shape will remain self-similar 
if the radial velocity v r at each point of the boundary is 
proportional to the radius at that point. The normal ve- 
locity v n — M*k is the projection of the radial boundary 
velocity onto the normal direction. The curvature is the 
derivative of the slope with respect to the arc length 



d9 _ 0' cos(0 - 9) 
ds r 2 d 



(Bl) 



Thus, the condition of the self-similarity of the shrinking 
shape can be written as 



M* = C r 2 d, (B2) 

r2d 



where C is some proportionality constant. To complete 
the description we need to express the radius r 2 d(0) in 
terms of i 



r' 2 M =r 2d (0) sin(0-0(0)). 



(B3) 



Without loss of generality we set r 2 d (0) = 1 and integrate 
equations (|B2|) and i|B3|> together up to (f> — ir/Q. The 
second boundary ^-boundary condition 8 (tt/6) = tt/6 se- 
lects a unique C. The numerical shooting yields C ~ 
0.903535. The shape of the self-similarly shrinking grain 
on a hexagonal lattice is remarkably close to a circle. The 
largest and smallest grain diameters differ by only 0.4%! 
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